Thermal spin transport and spin-orbit interaction in ferromagnetic/non-magnetic 

metals 



O 
(N 

1—5 



I 



a 



I 

o 
o 



> 

o 

rn 

i> 

O 



X 



Abraham SlachterQ Frank Lennart Bakker, and Bart Jan van Wees 
Physics of Nanodevices, Zernike Institute for Advanced Materials, University of Groningen, The Netherlands 

(Dated: July 19, 2011) 

In this article we extend the currently established diffusion theory of spin-dependent electrical 
conduction by including spin-dependent thermoelectricity and thermal transport. Using this theory, 
we propose new experiments aimed at demonstrating novel effects such as the spin-Peltier effect, 
the reciprocal of the recently demonstrated thermally driven spin injection, as well as the magnetic 
heat valve. We use finite-element methods to model specific devices in literature to demonstrate 
our theory. Spin-orbit effects such as anomalous-Hall, -Nernst, anisotropic magnetoresistance and 
spin-Hall are also included in this model. 



I. INTRODUCTION 

Spintronics uses the spin degree of freedom to demon- 
strate new functionality in ferromagnetic/non-magnetic 
hybrid devices^. In time, many new functional devices 
have been proposecP'S'and measured^ ^ utilizing the spe- 
cial properties of spin transport. Recently, the coupling 
between thermoelectricity and spin transport has been 
added to this field. New applications resulting from this 
coupling ar e su mmarized under the branch called spin- 
caloritronica^'SH^. 

A diffusive transport theory for spin-depe ndent elec- 
trical conduction is currently well establishecP^EIl. This 
theory has been extended to non-collinear systemsiS^ 
which becomes relevant when spin-dependent tunnel- 
ing through interfaces is considered or to quantify dy- 
namic processes such as spin-transfer torque^^' or spin- 
pumping'^''^. 

In this article, we extend the collinear theory of dif- 
fusive transport for spin-dependent conduction to in- 
clude spii^dependent thermoelectricit}-^^ ^'*, spin-orbit 
effect ^^'251261 g^j-^^ ^Yso spin-dependent thermal transport. 
We use finite-element methods to demonstrate our the- 
ory in order to extract useful para meter s from complex 
three-dimensional device geometriea^^Sl! Various recent 
experiments are taken from literature to extract the pa- 
rameters which govern the effects. 

The setup of this article is as follows. In section II, we 
begin with a description of finite-element modeling where 
we specify the structure of the model and the solvers 
used. In Section III we describe how to make finite- 
element-models which describe electrical spin-transport. 
We illustrate this model by calculating a recent exam- 
ple from literature's. We also show how the direct and 
inverse spin-Hall effect can be included and use it to 
model an experiment by Kimura et al.l^. A thermo- 
electric model which excludes spin is described in sec- 
tion IV wher e it is shown how spin-orbit effects can be 
includecPSEa. 

In section V, we introduce the thermoelectric-spin 
model. This model can describe the individual effects 
related to thermoelectricity or spin-dependent electrical 
transport. However, the introduction of spin-dependent 



thermoelectric coefficients also allows to demonstrate two 
new physical effects: the recently demonstrated thermal 



spin injection^ and its Onsager reciprocal effect: the 
spin-Peltier effect. Thermal spin injection describes the 
injection of spins in a non-magnetic material when a heat 
current is sent through a ferromagnetic/non-magnetic in- 
terface. The spin-Peltier effect describes spin-dependent 
heat transport across this interface due to the injection 
of spins in the ferromagnetic material. 

In section VI, a phenomenological theory for spin- 
dependent heat transport is proposed, where the con- 
cept of a spin temperature is introduced and the thermal 
analogy of the spin valve: the magnetic heat valve^. 
We apply the model on a previously measured samplers! 
to determine an upper limit for the relaxation of spin- 
dependent heat at room temperature. Thermoelectric- 
ity not only connects spin-dependent charge transport 
to heat transport but also connects spin-dependent heat 
transport to charge transport. This provides new ways 
to generate spin temperatures and to detect these. 

We conclude this article with a discussion on how 
the spin-dependent effects can be described by a the- 
ory which includes tunneling through interfaces or non- 
collinear magnetizationa^^l. 



II. FINITE-ELEMENT MODELING 

The finite-element modeling in this article is performed 
using the software package Comsol Multiphysics (version 
3.5). It solves partial differential equations (PDE's) for 1, 
2 or 3 dimensional geometries defined in a CAD drawing 
program. In a diffusive transport theory, the PDE's are 
determined by the conservation of the generalized cur- 
rents for the physics considered. These can be fo rmally 
derived from Boltzmann transport theorjEEil^ xhe 
fiuxes, put into a vector by J = (J„j, J„2, ...), are gov- 
erned by a vector of continuous variables u = (ui, U2, ■•■) 
through the conductance matrix c: 



J = ~cVu (1) 

Depending on the dimensionality n (ID, 2D or 3D) 



of the finite- element model, the elements of the fluxes 
Jui, Ju2i ■■■ are vectors themselves of size n. They deter- 
mine the currents in the respective directions defined by 
the coordinate system of the model. The elements of the 
conductance matrix C(i,j) are then n x n matrices. For 
an isotropic model, these are scalar matrices while for 
anisotropic transport the elements can be different. The 
PDE's in the bulk are determined by the conservation of 
fluxes: 



V • J = / (u) 



(2) 



Where a source term / (u) exists which may depend on 
the variables themselves. As an example, for simple elec- 
trical transport u = V , J = Jc, c = a and /=0. Here 
V is the voltage, J^ the charge current and a the electri- 
cal conductivity. Eq. [TJthen states Ohm's law while Eq. 
[2] is the Poisson equation representing the conservation 
of charge. The system under consideration is solved by 
stating the boundary conditions. These can be set for 
each variable (a Dirichlet condition) or flux (a Neumann 
condition) individually. In our example of electrical con- 
duction, a charge current can be sent through the mate- 
rial by setting the charge current to a specific value at 
one interface and the voltage to a specific value at an- 
other. The outer interfaces are insulating Jc—0 and the 
currents are continuous across internal interface Jc\i — 

Jch- 

A (tetrahedral) mesh of typical 300k elements is cre- 
ated by the finite-element program where specific de- 
tailed meshing is often used in the areas of interest by 
specifying a minimal element size. The PDE's are solved 
using a built-in (non-linear) solver which uses the itera- 
tive generalized minimal residual solver (FGMRES) with 
a geometric multigrid preconditioner, which in its turn 
uses a a direct sparse object-oriented linear equations 
solver (SPOOLES). 

The models we use are generally non-linear and the 
device on which the model is based is measured electri- 
cally. Therefore, the resulting measurable voltage is also 
non-linear: 




FIG. 1: (Color online) Results of the modeling of the non- 
local spin valve structure used by Yang^ et al. a) A spin 
current is injected into the Cu bar which connects two fer- 
romagnets FMi, FM2 by sending a charge current Ii_2 = 
1 mA. The injected spins diffuse towards FM2 where they 
are absorbed. The color shows a calculation of spin voltage 
Vs ^V^ — V^aX Ii_2 = 1 rnA for the structure with a 4 nm thin 
ferromagnet FM2. b) Calculated spin valve signal ¥3-4/11-2 
(mSl) versus the spin polarization for electrical conductance 
for a 20 nm thick FM2. The grey area shows the measured 
spread. 



III. SPIN TRANSPORT 

Spin-dependent electron transport in systems 
consisting of collinear magnetizations and clean 
ferromagnetic/non-magnetic interfaces is commonly 
described in terms of a 2-channel model. First suggested 
by MotifiS^ and later derived from the Boltzmann 
transport theorjff^^ it describes electrical conductance 
separately for spin-up (f) electrons, the component 
parallel to a magnetization, and spin-down electrons 
{D, the antiparallel component. Each channel has its 
own conductivity u^,^, voltage V^^^ and charge current 
Ji;.\- Usually a simplified resistor model is employed to 
describe spin-dependent transport. While it is sufficient 
for many approximations, it can become inaccurate for 
complex three-dimensional structuregSZI. 

Spin-dependent transport can be modeled by a set of 
PDE's by using the spin-dependent voltages as variables 
u = {V^,V\J). The spin dependent currents J = {J^,J\J) 
are then defined by the spin-dependent conductance ma- 
trix: 



V ^Ril + R2l^ + Rsl^ 



(3) 



The contributions i?„(y/A") to the nonlinear voltage 
can be separately determined from experiments. By us- 
ing multiple lock-in systems which are set to measure the 
different responses w, 2a;, 3uj resulting from a sinusoidal 
charge current / of frequency uj sent through the device, 
it is possible to determine these contributions^S. We may 
extract them from the constructed model by studying 
how the simulated voltages depend on the applied charge 
current. The nonlinear contributions RniV/A"') are de- 
termined by calculating the model at currents of ±7 and 
±2/ and deriving these contributions. Here / is typi- 
cally of a size used in experiments where all interesting 
contributions Rn{V/A"-) are considerable. 



a^ 
a^ 



(4) 



The conservation of charge current is given by V 
'jt + Ji) = 0- The Valet-Fert equation \/^{Vf - V^) = 



''"^2 '*' is derived from the conservation of spin currents. 
Defining a spin polarization for electrical conductance 
P/ = {a^ — (j^)/{a^ H- a^) we find a source term: 



/ = 



4A2 



(n - n) 



I 



(5) 



The inputs for this model are the specific geometry, 
the conductivity a — cr^ -I- cr^, spin relaxation length 



A and spin polarization Pj. These parameters can be 
determined from various experiments. For example, 
the relaxation length of non-magnetic materials can be 
determined by varying the distance between two fer- 
romagnet materials in a spin valve or measuring spin 
precessiorp3E2^ The relaxation length of fcrromagnets 
can be determined by angle-resolved photoemission^^, 
while the conductance polarization of fcrromagnets can 
be investigated by measuring the Doppler shift of spin- 
waveJ^. 

While the model typically works well for pillar struc- 
tures with clean interfacea^Slj in lateral systems the 
model generally overestimates the observable spin-valve 
signalaS^l. This is because it is often necessary to per- 
form ion-milling prior to deposition to obtain ohmic in- 
terfaces. This is taken into account in the modeling by 
using a reduced conductance polarization of the ferro- 
magnet P/. For example, for a permalloy (NigoFe2o) 
ferromagnet commonly used in a lateral spin valve, the 
polarization is reduced from 70 %, de termined from mea- 
sured Doppler shifts, to 30-509#'^. 

An example of the application of this model can be 
found in Bakker^. Another example of where this model 
can be applied is shown by_Yang221 et al. Here, the non- 
local spin valve geometrjI^SSH is used to inject a spin 
current Jj- — J^ free of charge current (a so-called pure 
spin current) into a thin ferromagnet. The magnetization 
of thi s ma gnet is switched by the resulting spin-transfer 
torqu^2120] ^ threshold exists for this process given in 
terms of the charge current which should be sent through 
the device. 

A model of this device geometry with resulting spin 
voltage is shown in Fig. [T] Using the measured conduc- 
tances, relaxation length of copper Xcu = l/irn and spin 
valve signals, we determine an effective spin polarization 
of P/ «0.6 from the measured spin valve signals 9-21 
TuQ from a batch of samples with thick fcrromagnets. 
By performing an integration of the spin current flowing 
through the FM2/NM interface, we find that the amount 
of pure spin current injected into the second ferromag- 
net Is is 13% of the total charge current / sent through 
the device. Using an effective formula for spin-transfer 
torque switching^S! we find that the ferromagnet should 
switch at Is=930/zA using common parametersSi where 
in the experiment the required charge current I=5mA 
results in a a spin current of Is=675/iA. Considering the 
empirical spread found for spin valve signals, this is very 
reasonable. 

Because the electrical current density spread through- 
out the device is modeled, the Biot-Savart law also allows 
to calculate the magnetic fields present in the device. The 
magnetic fields at position f is determined by performing 
a volume integral over the entire device: 



Biro) = ^ I ^'^'J'll'^ dr 



47r J \r — Tor 
Here /^o is the magnetic vacuum permeability. We find 



that the magnetic field at the center of the switchable 
ferromagnet is \B\ = 1.2 mT at the maximum applied 
charge current of 5 mA while the switching field of the 
magnet is 8 mT. This directly rules out that the magne- 
tization switches by the induced magnetic fields, but it 
can be responsible for observed asymmetries in the spin- 
transfer torque switching process. 



A. Spin-Hall Effect 

Spin-orbit effects in fcrromagnets are often sizeable due 
to their complex band alignment. As a result, measure- 
ments on spin-orbit effects in ferromagnets were already 
reported more then a century ago. First D'Yakanov and 
PereP^and later HirsclP suggested that the same process 
which governed these effects in ferromagnets, whether it 
be due to band alignment (intrinsic) or spin-dependent 
scattering (extrinsic), can also be responsible for a new 
effect in paramagnetic materials: the spin-Hall effectPSS. 

The direct spin-Hall effect describes that when a charge 
current Jc is sent through a material with strong spin- 
orbit interaction, a spin current Js flows away from the 
center of the conductor with its spin direction m, a unity 
vector, perpendicular to the charge and spin current. In 
the indirect version, a spin current flowing through the 
material creates a voltage perpendicular to the spin and 
spin current direction. Based on a Boltzmann transport 
theory derived by Zhang^il, the effects are governed by 
the following two equations: 



VV 



SH 



ysH 



m X Jg 



yylSH ^ ^o^^ ^ J 



a 
a 



(7) 
(8) 



Here \IV^^ and VV/'^^ are the bulk charge and spin 
voltages resulting from the direct and indirect spin-Hall 
effect and BgH is the spin-Hall angle, typically a small 
fraction of one. Both effects can be included into the 
2-channel model. To do this, we rewrite both equations 
into spin-up and spin-down currents and obtain the new 
spin-dependent conductance matrix: 



a 



SH „ 



(9) 



Where non-diagonal elements a^^ are included. These 
become skew-symmetric matrices determined by the 
spin-direction m considered in the device: 

c^t.f (*'■?') = TQsHO^\^ei3kmk (10) 

fe 

Here (i,j,k) are the indices of the predefined xyz axes and 
Eijk is the Cevi-Levita symbol. 
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FIG. 2: (Color online) Calculated spin voltages Vs (piV) 
for the device of Kimura^^ et al. a) In the direct spin-Hall 
configurations a charge current Ii_2 = 1 mA is sent from fer- 
romagnet to copper arm and the resulting spin-Hall voltage 
V3_4 is measured on the Pt strip, b) XZ cross section in the 
middle of the platinum strip in the indirect spin-Hall configu- 
ration. The charge current l3_4 = 50 /lA is short-circuited by 
the copper strip which is why most spin accumulation enters 
through the corners. When the spin voltage enters the cop- 
per strip it is only a small fraction of the ±8/xV spin voltage 
present in the bulk platinum. The voltage Vi_2 is now mea- 
sured by converting the spin voltage Va to a charge voltage 
at the FM/NM interface. 



To demonstrate the theory, we model a device mea- 
sured by Kimur£p3 et al. where both the direct and indi- 
rect spin-Hall effect were measured in a single nanoscale 
device at room temperature for the first time. The re- 
sults from this model are shown in Fig. [2] The device 
consists of a single permalloy ferromagnet which is con- 
nected to a 4 nm thin platinum strip by a copper cross. 
A pure spin current can be injected into the platinum 
strip by sending a charge current Ii_2 from the ferro- 
magnet to one of the arms of the copper cross. When 
the magnetization of the ferromagnet is aligned in the 
±y-direction, the spin current flowing into the platinum 
in the -z direction creates a voltage due to the indirect 
spin-Hall effect in the x-direction. This voltage can be 
measured between the two contacts present on the plat- 
inum strip. In the same device, sending a charge current 
through the platinum strip creates a spin current flowing 
in the z-direction where the spins are aligned in the ±y- 
direction, which is now due to the direct spin-Hall effect. 
When this spin current arrives at the permalloy strip, 
it is converted into a voltage which can be electrically 
measured. 

Both the direct and indirect spin-Hall signal in this de- 
vice is around 60 /lil at room temperature at a distance 
of 400 nm from platinum strip to ferromagnet. Using 
the common paramete rs Pp„= 0.3, Xpy = 5 nm, Apt — 2 
nm and Xcu = 350 nnpSEMZI ^nd the measured conduc- 
tivities, we find that a spin-Hall angle of 0sh = 5-10^^ 
accurately produces both signals. This angle is around 
8 times smaller then previously deduced. The difference 
lies somewhat in the used parameters, which are differ- 
ent, but mostly in the modeling which takes into account 
effects which the resistor model does not. For example. 



when a charge current is sent through the contacts con- 
nected to the platinum strip, the charge current mostly 
goes through the copper on top of the platinum strip 
instead of the platinum strip itself, which is due to the 
high conductivity of the copper. Copper is expected to 
have virtually no spin-Hall effect which reduced the ob- 
served voltage. We calculate that only 1% of the charge 
current goes fully through the platinum strip. A sim- 
ilar effect also occurs for the indirect spin-Hall effect. 
When a spin current is injected into the platinum strip, 
the produced voltage is also short-circuited by the cop- 
per strip. This example clearly demonstrates the need 
for a three-dimensional model to accurately predict the 
relevant parameters. 



IV. THERMOELECTRICITY 

Ferromagnetism of typical (metallic) ferromagnets 
used in spintronic experiments originates from the elec- 
tronic band structure. This band structure varies 
more strongly than those of the non-magnetic parts 
due to ferromagnetically induced shifts in the s- and 
d-bands. The Mott-formulgP^ describes that thermo- 
electricity is strongest when the electrical conductivity 
varies strongly around the Fermi energy. Since this vari- 
ation is dependent on the band structure, the ferromag- 
nets used in spintronic experiments are good thermoelec- 
tric materials^. A model of thermoelectricity is there- 
fore highly requi red to explain the voltages observed in 
nanoscale device^^SlSZi ^^^^ ^.^^^ g^jgQ ^g used to demon- 
strate new physics^. 

A bulk thermoelectric theory based o n Bol tzmann 
transport has been in use for half a centurjB^ESI i^ this 
theory, in addition to charge transport Jc, the heat trans- 
port Q [W/m^] is also taken into account. In a thermo- 
electric model, we use the variables u = {V, T) and the 
currents J = {Jc,Q)- The transport is then determined 
by the conductance matrix: 



a crS 
an k 



(11) 



where k is the thermal conductivity, S the Seebeck coeffi- 
cient and n the Peltier coefficient. The Peltier coefficient 
is related to the Seebeck coefficient by the thermody- 
namic Thomson-Kelvin relation H = ST, where T is the 
reference temperature. The Seebeck coefficient governs 
the Seebeck effect: the generation of a voltage as a re- 
sult of a temperature gradient. The Peltier coefficient 
describes the transport of energy when a charge current 
is sent through a material. 

In the model, the charge current remains conserved. 
However, the heat current is not. When a charge cur- 
rent is sent through a material, the inelastic scattering of 
electrons causes a rise in temperature. This is a process 
known as Joule heating. Conservation of the currents 
gives the source term: / = (0, J^/cr). 



Unlike electrical transport, which is only carried by 
(missing) electrons, thermal conductivity typically has 
two contributions: that from the energy transported by 
electrons and from lattice vibrations: phonons. In met- 
als, the transport by electrons often dominates. In this 
case, the Wiedemann-Franz law k = aLT, which relates 
thermal and electrical conductivity through the Lorenz 
number L, can be used to estimate thermal conductivi- 
ties from measured electrical conductivities. 

This model has recen tly b een demonstrated in 
nanoscale spintronic device^^filSZI, These experiments il- 
lustrate that the ferromagnetic/non-magnetic junction 
used acts as a thermocouple to measure the local temper- 
ature of the device. In addition, when a charge current 
is sent through such a junction it heats or refrigerates 
the device due to the Peltier effect. In the non-local spin 
valve geometry, the interplay between both effects leads 
to a baseline resistance i?i which is purely thermoelectric 
in nature. 

The inclusion of Joule heating makes heat transport, 
and therefore also charge transport, nonlinear. This re- 
sults in nonlinear voltage signals. While the linear re- 
sponse Ri is often resulting from conventional Ohmic 
path^^, the higher harmonic response are the result of 
thermoelectric effects. For example, the i?2 response 
can result from Joule heating of devices measured by 
a ferromagnet/non-magnetic thermocouple. If the typ- 
ical bulk thermal dependencies 7, a of the Seebeck co- 
efficient S(T)=S'o(l -f- 7Ar) and conductivity a{T) = 
cro/(l -I- a AT) are used in the model, it is possible 
to explain even higher order nonlincaritics observed in 



experimeni 



Ea 



A. Spin-Orbit Effects 

Whenever magnetic fields are applied or magnetism is 
involved in experiments, a large variety of spin-orbit ef- 
fects may occur. These effects break the symmetry of 
the thermoelectric model. The most notable of these is 
the Hall effect which describes the generation of a volt- 
age perpendicular to the plane made by an applied charge 
current and magnetic field. The observed transverse volt- 
ages generated in a ferromagnet are typically larger than 
expected from the ordinary Hall effect resulting from in- 
ternal magnetic fields. Therefore, in ferromagnets, it is 
named the anomalous-Hall effect. The exact origin of 
this effect has been under continuous debate for over a 
century and seems to depend on the specific conditions 
of the experimenlPZl. 

The equation which describes the anomalous-Hall ef- 
fect is obtained from Boltzmann transport theorySl 
^yAH _ _0^^ff^ X j^ where Jah is the contribution 
to the current due to the anomalous-Hall effect. It is 
incorporated in the model using a skew symmetric con- 
ductivity: 



o'(ij) 



«■ I Sij - 9 AH 2J ^ijkmk I 



(12) 



Where 6ah is the anomalous-Hall angle. Such contribu- 
tions are often called spurious in spintronic devices since 
they can mimic spin valve signals ■^■^. However they can 
be modeled in spintronic devices to calculate their mag- 
nitude. 

The Seebeck, Peltier and thermal conductivity coeffi- 
cients have similar spin-orbit contributions. The effects 
which arise due to these contributions are named ther- 
momagnetic effects. The effect arising from the Seebeck 
coefficient is named the anomalous-Nernst effect^ and 
describes the generation of a voltage perpendicular to 
the plane made by an applied heat current and magne- 
tization. The contribution to the Peltier effect is named 
the anomalous-Ettinghausen effect^^ while the contribu- 
tion to the thermal conductivity is named the anomalous 
Righi-Leduc effecfP. 

In addition, spin-orbit effects also determine that the 
conductivity of a ferromagnet is different by measuring 
it parallel or perpendicular to the magnetization of the 
ferromagnet. This is the anisotropic magnetoresistance 
effect and can be included in the model by adding an 
anisotropic contribution to the conductivity. The sym- 
metric contribution then becomea^Sl 



^■(U) = 0-_L (% - RAMRmiTTlj) 



(13) 



Where 5ij is the Kronecker delta and Kamr the co- 
efficient governing anisotropic magnetoresistance, typi- 
cally a few percenlP. By symmetry, we expect equiva- 
lent anisotropic contributions to the Seebeck, Peltier and 
thermal conductivity. However, specific measurements 
demonstrating these effects have not been reported to 
the authors knowledge. 

We recently measured two of the previously discussed 
effects in a spin-caloritronic device^Sl. This device is best 
described as the thermal equivalent of a spin- valve struc- 
ture (the magnetic heat valve) and will be discussed 
in more detail afterwards. Here, one of the ferromag- 
nets was Joule heated to generate a heat current Q 
through the spin-valve device. The temperature was 
measured on the second magnet using a thermocouple. 
Anisotropic magnetoresistive heating and anomalous- 
Nernst were found to dominate the measured voltage 
behavior. By virtue of a finite-element model, it was 
possible to determine magnetization angles and also the 
size of the anomalous-Nernst effect. 



V. THERMOELECTRICITY AND SPIN 

Thermoelectricity extends charge transport theory and 
includes effects governed by the Seebeck, Peltier and 
thermal conductivity coefficients. The spin-transport 



model extends charge transport theory to include the 
spin-dependency of the conductivity and introduces 



the concept of spin-dependent voltages V 



T,;- 



The 



model which extends charge transport theory to include 
both thermoelectricity and spin-transport is named the 
thermoelectric-spin model. It has been used for almost 
50 years to describe thermoelectricity in ferromagnets^^ 
and more recently, to describ e thermoelectricity of 

multilayered spin valve d^^ l '^^ l ^^ l and spin transport in 
ferromagnetplHna] 

The relevant physics and measurable voltages in de- 
vices can be calculated using finite-element modeling. 
The spin dependent voltages and temperature are the 
variables u — (V^,y^,T) and the fluxes are determined 
by the spin-dependent charge currents and heat current 
J = (J^, J^,Q). The conductance matrix now allows us 
to include spin-dependent Seebeck S-t-,^ and Peltier co- 
efficients n^.^ to describe not only the coupling between 
charge and heat transport, but also the coupling between 
spin and heat transport. The conductance matrix is given 
by: 



(7^ (T^ S*-^ 



(14) 



The conservation of charge currents remains un- 
changed. However, the Valet-Fert equation is altered be- 
cause in the derivation thermoelectricity is disregarded. 
It is originally derived using particle conservatioiPSl 



e 



t,i 






(15) 



Here t^^ represents the time for a spin up electron to flip 
its spin to spin down while t^^ represents the time from 
a spin down electron to flip its spin to spin up. The ex- 
cess electron densities are given by the Einstein relation 
^ti — N^i^y^i^ with N^^ the spin dependent densities 
of states at the Fermi energy. In the thermoelectric- 
spin model, the spin-dependent charge currents J^^^ = 
— cr-|-^^ [W + S^^^VT) additionally includes a tempera- 
ture gradient as well as the spin-dependent Seebeck co- 
efficients. 

The Seebeck coefficient describes how the conductivity 
depends on energy and is described by the Mott formula. 
By virtue of the Einstein relation, the Seebeck coefficient 
is determined by the energy derivative of the density of 
states (^)p and the relaxation time (jg)p at the 
Fermi energy. To develop an altered Valet-Fert equation, 
the energy dependence of the densities of states N^^^ and 
relaxation times t^\^.\.-\ needs to be taken into account 
at the right side of Eq. [T5j While theoretically these 
contributions can be taken into account, in practice not 
much is known about these specific energy dependencies. 
For simplicity, we ignore such terms in the modeling and 
note that they can be responsible for small bulk source 



term9^BIl_ Conservation of charge, spin and heat currents 
are now taken directly from the individual thermoelec- 
tric and spin-transport models and produce the following 
source term: 



/ = 



(1-P>; 









(16) 



Where we introduced the individual Joule heating of both 
spin channels J^ ,/a^^i. 

In the spin-dependent charge transport model, the 
ferromagnetic/non-magnetic interface plays a crucial role 
in converting spin transport into charge transport and 
vice versa. Using this coupling between spin and charge, 
magnetic memory elements can be constructed. In 
the thermoelectric-spin model, the ferromagnetic/non- 
magnetic interface plays a similar role. In the following 
we will show that this interface can be used to convert 
heat transport into spin transport and also vice versa. 
The conversion of a heat current into a spin voltage de- 
pends solely on the spin-dependency of the Seebeck coef- 
ficient and is therefore named the spin-dependent See- 
beck effect or by its more application oriented name: 
thermal spin injection. This was recently measured in 
a dedicated spin-caloritronic device23. 

Here, we also propose a measurement scheme for the 
reciprocal effect. We show that when a spin current is 
injected into a ferromagnet a net heat fiow develops, even 
in the absence of a charge current. This effect, named the 
spin- Peltier effect, depends solely on the spin-dependency 
of the Peltier coefficient. 



A. Spin-dependent Seebeck effect 

The spin-dependent current model dictates that when 
a charge current J^ = J^ + J^ is sent through the bulk of 
a ferromagnet, a spin current Js = J^ — J^ accompanies 
it of which the size is determined by the conductivity 
polarization Js ~ Pi Jo- 
in similar fashion, the thermoelectric-spin model dic- 
tates that when a heat current Q is sent through the 
bulk of a ferromagnet in the absence of a charge current, 
a spin current Js = —api^ — P^)S'sVT/2 flows, of which 
the size is determined by the spin-dependent Seebeck co- 
efficient Ss = S^ — S^ = PgS. Here Ps is a fraction of 
the regular Seebeck coefficient S, defined in terms of the 
spin-dependent Seebeck coefficients as: 
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(17) 



In a non-magnetic material, both the conductivity 
polarization and the spin-dependent Seebeck coefficient 
are zero. When a charge current is sent through a 
ferromagnetic/non-magnetic interface, the discontinuity 




FIG. 3: (Color online) Thermal spin injection, measurement 
scheme and measured resullf^. a) A heat current sent through 
a F/N interface generates a spin voltage 14 ~ XfSsVT at 
the interface which extends a distance Xm, Xf in the mate- 
rials, b) Thermal spin injection can be measured by Joule 
heating FMi in a lateral spin valve. This generates a heat 
current Q over the FMi/NM interface which injects spins. 
The spin voltage is turned into a measurable voltage by the 
FM2/NM interface. The size of thermal spin injection can be 
determined by selectively switching the magnetizations, c) 
Measurement result. A signal due to thermal spin injection 
is present on a large background caused by the measurement 
of Joule heating by the FM2/NM thermocouple. In addition, 
small traces of anisotropic magnetoresistive heating effects as 
well as anomalous-Nernst effects are present. These can be 
seen by the small dip at the switching field of FMi and the off- 
set in background voltage between large positive and negative 
magnetic fields (see also fig. 7). 



interface is measured by a second ferromagnetic/non- 
magnetic interface which converts the spin voltage into a 
measurable voltage. A schematic picture of thermal spin 
injection and the used measurement scheme is shown in 
Fig. [31 

Because Joule heating scales quadratically with the ap- 
plied charge current /^, thermal spin injection results in 
a nonlinear spin dependent signal R2 — R2 — R2 where 
R2 and i?^^ are the parallel and antiparallel contribu- 
tions. The measured result is shown in Fig. [3J;. The 
applied temperature gradient was very limited due to the 
relatively large lateral size and because electromigration 
prohibits heating in excess of 40K in this particular case. 

In addition to thermal spin injection, we observe small 
traces of spin-orbit effects such as the anomalous-Nernst 
effect in the second ferromagnet and anisotropic mag- 
netoresistive heating of the first ferromagnet. These ef- 
fects have been more thoroughly examined in another 
device^SI. These effects express themselves by a differ- 
ence in background voltage for both parallel orientations 
and the observed small curve in voltage prior to the low 
field switch. 

An application of spin currents lies in its ability to 
switch the (uniform) magnetization of a ferromagnet 
around its easy axis by means of spin-transfer torque^^l. 
This effect has a threshold in the spin current which needs 
to be injected in a small volume ferromagnet. Since elec- 
trical and thermal spin injection have the same physical 
origin, the discontinuity in bulk spin current, we may di- 
rectly compare the critical temperature gradient in the 
ferromagnet needed to switch a F/N/F spin valve by spin 
transfer torque to the critical charge current density for 
this switching process. If the critical charge current den- 
sity is known, we can calculate the critical temperature 
gradient which is needed for the switching process: 



in bulk spin current creates a spin voltage at the inter- 
face and injects a net spin current in the non-magnetic 
material. The same situation occurs when a heat current 
in the absence of a charge current is sent through the 
interface. 

This effect, the injection of spins in a non-magnetic 
material by a heat current, is named thermal spin injec- 
tion or, equivalently, the spin-dependent Seebeck effecP^. 
Since both electrical and thermal spin injection arise 
from the discontinuity of the bulk spin currents, both 
effects have similar behavior. For example, the spin volt- 
age spreads an equal distance away from the interface 
and both effects suffer from the conductivity mismatch 
probleirPSI which strongly reduces spin injection in low 
conductivity materials such as semiconductors. 

Thermal spin injection was recently demonstrated in 
a multiterminal lateral device^S. In this device, a tem- 
perature gradient is applied to a F/N/F spin valve by 
Joule heating one of the ferromagnets with the help of 
a large charge current. The thermoelectrically generated 
spin voltage across the first ferromagnetic/non-magnetic 
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Here Jcxit is the threshold in charge current at which 
spin transfer torque switching takes place, a the ferro- 
magnetic conductivity and VTjcrit the critical tempera- 
ture gradient in the bulk ferromagnet. As an example, if 
we assume a critical charge current density Jcrit — 10^^ 
A/m^ for a permalloy ferromagnet, with the common (es- 
timated) parameters P/ — 0.6, Ps — 0.6, a — A- 10® S/m 
and S ~ — 20 fiV we find a critical temperature gradient 
of VT = 4 • 10^ K/m. A typical F/N/F stack of 25 nm 
then switches at an applied temperature difference of 100 
degrees. This process is known as thermal spin transfer 
torque, and recently evidence has been found for it^^. 

The ability to use finite element modeling should al- 
low to engineer multiterminal F/N/F pillar devices which 
switch by thermal spin transfer torque. Such devices can 
combine the high polarization properties of pillar devices 
with the flexibility of lateral devices. By selectively heat- 
ing the device, the effect can also be used to lower the 
effective threshold in electrical spin injection. 
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FIG. 4: (Color online) Schematic representation of the Spin- 
Peltier effect. A spin current Js — Jf — J], free of charge 
current (J-|- + Jj. = 0) is injected from the non-magnetic side 
of the F/N junction into a ferromagnet FM. The top shows 
the resulting spin-dependent voltages, calculated using the 2- 
channel model. Despite the fact that no charge current is 
flowing through the junction, a net heat flow Qn ~ |(nt — 
Dj.) Js develops in the ferromagnet which quickly drops off due 
to the spin relaxation length \p. Depending on the sign of 
the spin current and the parallel/antiparallel alignment of the 
magnetization, net heat is transported from the non-magnetic 
material to the ferromagnet or vice versa. This creates a 
temperature difference AT between the bulk non-magnetic 
material and the bulk ferromagnet. 



B. Spin-Peltier effect 

The Onsager reciprocity relation dictates that when 
heat transport induces spin transport free of charge cur- 
rents in a ferromagnet, the opposite can also occur. A 
pure spin current injected in a ferromagnet should induce 
net heat transport. This reciprocal effect is named the 
spin-Peltier effect. 

We illustrate this effect in Fig. |4] by considering the 
F/N interface previously used. A pure spin current is 
injected from the non-magnetic side into the ferromag- 
net. When the spin current enters the ferromagnet it 
reduces in size at the spin relaxation length, for metals 
ranging from a few to tens of nanometers. The spin- 
dependent voltages which result from the spin current 
are sketched in the top part of Fig. [4J In the absence 
of charge currents, the heat current in the system due 
to the spin-Peltier effect Qn can be deduced from the 
thermoelectric-spin model: 



?n = ^(nt-n^)Js 



(19) 
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In the non-magnetic material, 11^ = 11^ and no net 
heat transport due to the spin-Peltier effect takes place. 



FIG. 5: (Color online) Proposed non-local measurement 
scheme for the spin-Peltier effect, a) The non-local spin valve 
geometry allows to inject a pure spin current into ferromag- 
net FM2 by sending a charge current Jc across the FMi/NM 
interface. This spin current creates a temperature difference 
across the FM2/NM interface changing the local tempera- 
ture Tp of the ferromagnet. The effect is detected by con- 
verting the temperature to a voltage Vtc = ASfa/Tf us- 
ing a thermocouple. Here ASfm is the difference between 
Seebeck coefficients between the ferromagnet and the non- 
magnetic material (yellow). The thermocouple measures the 
temperature Tp of the ferromagnet and optionally, makes use 
of the large Seebeck coefficient of the ferromagnet itself, b) 
The simulated spin-Peltier signal. The resulting signal should 
have a background determined by Peltier cooling/heating at 
the FMi/NM interface and subsequent consequence to the 
temperature Tp. A small signal of 100 nV due to the spin- 
Peltier effect should arise which only depends on the parallel 
or antiparallel alignment of both magnetizations. 



In the ferromagnet, the spin-Peltier coefficient, defined 
as lis = n^ — n^ can be non-zero. Close to the interface, 
the spin current Jg is also non-zero. This induces heat 
transport due to the spin-Peltier effect which drops off 
in the ferromagnet at the spin relaxation length. As a 
result, a temperature difference AT develops between the 
bulk non- magnetic material and the ferromagnet. 

The Thomson-Kelvin relation 11 = ST relates the con- 
ventional Seebeck and Peltier coefficients. This also holds 
for the individual spin species. From the recently found 
spin-dependent Seebeck coefficient Sg we can calculate 
the spin- Peltier coefficient Hg = SgT which can be used 
to estimate the effect. 

We calculate the exact temperature difference by con- 
sidering the total heat current Q = HsJs — kVT in the 
ferromagnetic and non-magnetic regions. Like spin and 
charge currents, the heat current is continuous across 
the interface. If we ignore Joule heating due to spin 
currents, it is also continuous in the bulk of the non- 
magnetic and ferromagnetic parts and equal to ki\/T\i. 
Here the index i denotes the region. The presence of 
the spin-Peltier effect induces an additional temperature 
gradient VT|n(a;) = j^Js{x) in the small ferromag- 
net region in which a sizeable spin current exists. The 
spin current drops off exponentially from the interface: 
Jsix) = JgC"^/'*'^. Here Jg is the spin current at the in- 
terface. If we integrate this additional temperature gra- 
dient over this region, we find the temperature difference 
between the bulk ferromagnet and non-magnetic mate- 
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This temperature difference depends solely on the spin- 
Peltier coefficient XI^, the spin relaxation length Xp and 
the thermal conductivity kp of the ferromagnet. Its sign 
is determined by the sign of spin current and the spin- 
Peltier coefficient. 

The non-local spin valve geometry is an ideal geometry 
to inject pure spin currents into a ferromagnet. The gen- 
erated temperature difference over the interface can be 
detected by measuring the temperature of the ferromag- 
net in which the pure spin current is injected. This can 
be achieved by placing a thermocouple on the ferromag- 
net. This measurement geometry is illustrated in Fig. [5] 
The background voltage is then solely determined by the 
Peltier heating of the FMi/NM interface which injects 
the spin current and the subsequent measurement of the 
temperature by the thermocouple. The spin-Peltier sig- 
nal then appears as a regular resistance Ri which solely 
depends on the parallel or antiparallel alignment of both 
magnetizations. 

For permalloy, all parameters are known^ and we can 
estimate the temperature difference which can be created 
in this manner. At a realistic maximum pure spin cur- 
rent which can be injected into a permalloy ferromagnet 
in the non-local spin valve geometry (see the previous 
discussion below Fig. [l]) we have J° — 10^^ A/m^. Us- 
ing this value, we find a temperature difference of AT 
= 20 mK between the parallel and antiparallel orienta- 
tion of the permalloy spin valve across the interface of 
the ferromagnet. A typical thermocouple which can be 
realized on a (lateral) ferromagnetP^ has an efficiency of 
AS* = AQfiV/K. This results in a maximal spin-Peltier 
signal of 800 nV which is small but observable. Initial ex- 
periments show signs of the spin-Peltier effect, however, 
it is hard to distinguish it from small parasitic effects, 
for example, the pick-up of regular non-local spin valve 
signals by an uneven distribution of the spin-voltage at 
the detecting interface. 



VI. BEYOND THERMOELECTRICITY AND 
SPIN: THE SPIN-DEPENDENT HEAT MODEL 

In the thermoelectric-spin model a single electron 
temperature was introduced which holds for both spin 
species. The energy of the electrons is distributed 
in their respective bands according to a position- 
dependent Fermi-Dirac distribution /|-.^(e, n^,r) with a 
spin-specific local density n^^\^ and local temperature T. 
This model assumes strong inelastic interaction within 
each spin species and also between them to obtain the 
required thermodynamic distribution. It is caused by 
electron-electron interaction or mediated by phonons 
through electron-phonon interaction. 



At low temperatures inelastic scattering becomes 
weaker and this requirement does not hold. It was shown 
in the past that inelastic scattering can be weak on 
the scale of the spin relaxation length in non-magnetic 
metals^ at sub-4K temperatures where electron trans- 
port is still diffusive, limited by elastic scattering. Al- 
though in this situation it is hard to speak of electrons 
which are distributed according to a Fermi-Dirac distri- 
bution in their respective bands, it is still possible to 
describe thermal transport according to a diffusion equa- 
tion. The temperature T then represents the local aver- 
age excess energy of electrons compared to the situation 
at zero Kelvin. In addition, the spin-dependent electron 
species also do not exchange energy with each other. 

This requires the introduction of a spin-dependent 
heat model where both spin channels have their own 
heat current Q^.^, thermal conductivity k^^\^ and spin- 
dependent temperature T^,\- This opens up the possibil- 
ity to demonstrate new thermal and, possibly, thermo- 
electric experiments in magnetoelectronic devices. Such 
a model has been first described by Heikkila et al^^. We 
introduce a bulk diffusion model and calculate an exam- 
ple device. 

The Wiedemann-Franz law describes the relation be- 
tween charge and thermal conductivity. Both are dom- 
inated by electron transport in metals. Therefore, a 
spin polarization in the electrical conductance Pi leads 
to a spin polarization in thermal conductance Pq, sim- 
ilarly defined in terms of spin-dependent heat conduc- 
tances as Pq = {k^ — k^)/{k^ + k^). The model for 
spin-dependent electrical and thermal transport now in- 
cludes the spin-dependent voltages and temperatures 
u = (V^, V^,T^,T^). The spin-dependent charge and heat 
currents J = [J^, J^, Q-^, Q^) are determined through the 
4x4 conductance matrix: 



(21) 



Where the spin-dependent thermoelectric effects, rep- 
resented by the coefficients S^^i and n^_^, are used in the 
relevant spin-dependent currents. The conservation of 
spin and charge currents remains the same, and therefore 
the components of the source term for the spin-dependent 
charge currents as well. It is straightforward to include 
the conservation of the total heat current Q = Q^ + Qi 
into its spin-dependent parts: the Joule heating of each 
channel J^ ,/aI , simply applies to the channels individ- 
ually. We note here that strictly speaking, if there is no 
inelastic scattering, there is no Joule heating. However, 
any weak inelastic scattering does raise the average en- 
ergy of the electron baths which allows Joule heating to 
be used in the model as a local source of heat. 

Although it is beyond the scope of this article to derive 
the conservation of spin heat currents Qs = Q^ — Qi 
from Boltzmann transport theoryl^l^ -^g may introduce a 
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phenomenological relaxation analogue to the relaxation 
of the amount of spins themselves, represented by the 
Valet-Fert equation for spin voltage. 

The difference in excess energy between both spin 
species is represented by the spin temperature Tg — 
T^ — T^. In our model, we assume a thermal equiva- 



lent of the Valet-Fert equation V T^ 



A? 



Here \q 



is the relaxation length for the spin temperature. This 
relaxation length is not only limited by spin flip pro- 
cesses, but can also be limited due to inelastic scatter- 
ing between both spin species, where energy is being ex- 
changed between both spin species without flipping its 
spin. This results in the boundary condition for spin re- 
laxation lengths \q < A. The spin relaxation lengths are 
equal whenever inelastic scattering is absent. 

New source terms ± — j^ — (T^ — T^) are then added 

such that the conservation of spin heat is also included. 
This leads to the following source term in this model: 
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The thermoelectric coefficients in this model typically 
scale with temperature and are very small at the tem- 
peratures where this model is applicable. For exam- 
ple, for many non-magnetic metals the Seebeck coeffi- 
cient scales linearly with temperature such that typical 
Seebeck coefficients are in the order of 10 nV/K at He- 
lium temperatures^. This also holds for typical ferro- 
magnets such as Cobalt or permalloy. Therefore, we first 
explore the special properties of spin-dependent heat con- 
duction itself and propose the thermal equivalent of the 
spin valve. 



A. Magnetic heat valve 
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FIG. 6: (Color online) The spin-dependent heat model, a) 
Equivalency between the coefficients of the spin-dependent 
heat and charge models when thermoelectricity and Joule 
heating is disregarded, b) A heat current Q sent through the 
F-N interface creates a spin-temperature Ts and spin-related 
temperature difference AT. c) The F/N/F magnetic heat 
valve. A heat current sent through a F/N/F spin valve struc- 
ture creates a temperature difference across it dependent on 
the parallel or antiparallel alignment of both magnetizations, 
d) A possible experimental realization of the F/N/F heat 
valve. By heat sinking one ferromagnet and Joule heating an- 
other, heat can be transported through the heat valve. The 
temperature of the second ferromagnet depends on the spe- 
cific thermal resistance of the heat valve, determined through 
the parallel or antiparallel alignment of the magnetizations, 
and can be measured using a thermocouple. 



which relaxes in the materials with the spin heat relax- 
ation length Xq. This is depicted in Fig. l6b. The size 
of the spin temperature at the interface can be directly 
deduced from the equivalently calculated spin voltage Vg 
in the charge transport modeP^I 



^ _ '^PqRf,qRn,q 
Q ^ Rf.q + (1 - Ph)RN,Q 
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If we disregard thermoelectricity as well as charge cur- 
rents in magnetoelectronic devices, such that Joule heat- 
ing is absent, spin-dependent charge and heat transport 
are each represented by an independent set of equations. 
The mathematical model introduced to describe spin- 
dependent heat transport is then identical to that which 
describes spin-dependent charge transport. The differ- 
ence between the models is the size of the coefficients. 
The equivalence of the coefficients used in both models 
is depicted in Fig. |6^. 

Consequently, concepts which are relevant in the spin- 
dependent charge transport model will have their equiva- 
lent in the spin-dependent heat transport model. A simi- 
lar resistor model also applies^ . For example, consider a 
heat current Q sent through a F/N interface. This creates 
a difference in temperature between both spin species Ts 



Where Rn,q = -^r^ and Rf,q — -f^ are the equivalent 
thermal resistances determined by the spin heat relax- 
ation lengths and the thermal conductivities of the ma- 
terials. A spin related 'thermal resistance' AT = ^PqTs 
also develops across the interface. 

There also exists a thermal equivalent of the electrical 
F/N/F spin valve. When a heat current Q is sent through 
a F/N/F spin valve, a temperature difference AT devel- 
ops across it, which depends on the parallel or antiparallel 
alignment of the magnetizations. We refer to this concept 
as the magnetic heat valve. It is depicted in Fig. [6];. In 
the electrical spin valve, a simple calculatiorPS gives the 
difference between parallel and antiparallel resistance per 
unit area Rp - Rap = 2P]RfRn/{Rf + Rn{1 - Pf)) 
whenever the distance L between both ferromagnets is 
L <C Xn- Whenever L ^ ^n,q we obtain the tern- 
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perature difference between the parallel and antiparallel 
alignment in the magnetic heat valve: 



ATp - ATap 

Q 



2PqRf.qRn,q 
Rf,q+Rn,q{1-PX, 



(24) 



As an example, let us consider a Py/Cu/Py heat valve 
in a 25 nm thick pillar stack where the non-magnetic 
metal is thin enough to satisfy the condition L <C Xn,q- 
Whenever the spin valve is held at a total tempera- 
ture where inelastic scattering is small but not negli- 
gible, say Xq ~ A/2 and assume the estimated values 
Pq == 0.6, Af = 5 nm, Ajv = l/xm, kcu = 300 W/m/K, 
kpy = 30 W/m/K an applied temperature gradient of 
lOK / 25 nm produces a significant temperature differ- 
ence of «2K across the spin valve depending on the par- 
allel or antiparallel alignment. 

In a spin valve, it is possible to use Joule or laser heat- 
ing to produce heat currents through a device^. How- 
ever, measuring a temperature difference across a device 
is non-trivial. Nevertheless, by measuring the absolute 
temperature at the second ferromagnet with the aid of a 
thermocouple, the process can be measured since the pro- 
cess does influence the net heat flow through the device. 
This experimental measurement technique is sketched in 
Fig. [BH. It is fairly non-trivial to use analytical solutions, 
as heat transport through a substrate and Joule heating 
itself are hard to calculate in a three dimensional geome- 
try. However, by fitting the obtained measured voltages 
to those resulting from a finite-element model with vary- 
ing geometry, it should be possible to extract useful coef- 
ficients, such as spin heat relaxation lengths at different 
temperatures. 

This experimental measurement technique has been 
used previously at room temperature^®. Here, a 
Py/Cu/Py spin valve was used and a Py/NiCr thermo- 
couple. A SEM picture of the device is shown with a 
typical measurement in Fig. [Tj The system was mod- 
eled with a regular thermoelectric model which showed 
that a maximum heat current of Q « 10° W/m^ can be 
achieved at which the temperature difference across the 
spin valve is w2K with a used charge current of /c = 2 
mA. 

Although at room temperature, regular spin-orbit ef- 
fects such as the anomalous-Nernst (1) and anisotropic 
magnetoresistance (2) are dominant, we can use this sam- 
ple to demonstrate a calculation of the spin-dependent 
heat model. Using the additional parameters shown in 
Fig. ItI and assuming no inelastic scattering {Xq = A) 
and pQ = P/, we calculate a temperature difference 
ATp — ATap — 8 mK due to the heat valve effect which 
leads to a response R| = 19.6 nV/mA^ when the ther- 
mocouple is measured. The temperature difference is 10 
times lower than the value calculated from Eq. [24] with 
VTIf ~ 3 • 10^ K/m and is due to spin relaxation in the 
non-magnetic material. With a noise level of sa5 nV/mA^ 
the heat valve effect should be observable if inelastic scat- 
tering between both spin species is absent. The absence 
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FIG. 7: (Color online) Modeling of a fabricated device^ 
which potentially could show the magnetic heat valve effect, 
a) SEM figure and measurement geometry. Two permal- 
loy ferromagnets (blue) are connected by a copper rectangle 
(brown) . The temperature Tf of the second ferromagnet is 
measured by a NiCr (contact 4,6) - Py thermocouple, b) Mea- 
surement at a typical current of 1 mA. No regular spin valve 
signal is present, which would be the result of the magnetic 
heat valve effect. The anomalous-Nernst (1) and anisotropic 
magnetoresistive heating (2) effects are present, c) Calculated 
spin temperature at I1-2 ~ 1 rnA zoomed at the copper rect- 
angle calculated using the parameters Xpy^q — 5nm, Xcu,q 
= 350 nm and Pq = 0.25. 



of a spin heat signal above the noise level shows that in- 
elastic scattering is strong enough such that we find the 
requirement Xq < 2^ valid at room temperature. When 
Aq = ^A the calculated signal is approximately identical 
to the noise level. 



B. Thermoelectricity and spin-dependent heat 

The spin-dependent heat model becomes relevant when 
inelastic interaction between the spin species is weak, 
which occurs at low temperatures. Thermoelectric ef- 
fects are small at these temperatures, which is why thus 
far we did not consider the connection between the spin- 
dependent thermoelectric effects and the effects due to 
spin-dependent heat. 

However, prospects in fabrication which connects the 
fiexibility of a multiterminal lateral device design with 
the high signals observed in pillar structures and the low 
noise experiments associated with a low operating tem- 
perature should increase the observability of the effects 



12 



so far considered. We may then also consider the higher 
order effects related to this connection. 

Whenever a charge current J^ is sent through a fer- 
romagnet a spin heat current Qg =^ HgJc also flows, de- 
termined by the spin-Peltier coefficient lis. Here we as- 
sumed Vf — Vi- Similar to the case of electrical spin 
injection and thermal spin injection, when this charge 
current is sent through a F/N interface, this creates a 
spin temperature T^ at the interface which relaxes in 
the respective materials at the respective spin heat re- 
laxation lengths Xq. We propose to name this effect the 
thermal spin-Peltier effect. The size of the effect is given 



by Eq. 23 with the source of the spin heat current in the 
bulk Qs = PqQ substituted by that due to this effect 
Qs = lis Jc- We note that this ignores the generation of 
a spin voltage at the interface which by ordinary thermo- 
electric effects is converted to a spin temperature. 

Although at low temperatures lis can be very small, 
the maximum charge current, for small devices typically 
limited by electromigration ( J™"'^ = 10^^ A/m^), is often 
larger than the typical heat current, which might render 
this effect more efficient to generate a spin temperature 
in a non-magnetic material than the previously described 
effect in the magnetic heat valveSS. 

The Onsager reciprocal effect can also occur in this 
model. Whenever a spin heat current Qs is injected into 
a ferromagnet at a F/N interface this creates a voltage 
difference AV\f~n between the bulk non-magnetic ma- 
terial and ferromagnet. We propose to name this the 
thermal spin-dependent Seebeck effect. The calculation 
of this voltage goes similar to the calculation of the tem- 
perature difference AT in the spin-Peltier effect. Check- 
ing the symmetry between the spin-dependent charge and 
heat models we can directly substitute the various coeffi- 
cients in Eq. [20] to obtain the induced voltage difference 
over the interface: 



AV^I 



F-N 



2kF 



Xf,qQs 



(25) 



It is possible to measure this voltage difference directly 
over an interface in a multiterminal non-local device. 
However, it requires a source of pure spin heat current 
Qs in which preferably a charge-related spin current J^ is 
absent. This is a situation difficult to achieve. However, 
if both sources of Qs and Js scale differently with ap- 
plied charge current it is perhaps possible to distinguish 
between the generation of a voltage over a F/N inter- 
face due to charge-releated spin currents and those due 
to spin heat currents in a suitably designed experiment. 



VII. DISCUSSION 

Throughout this article we have considered transpar- 
ent (Ohmic) interfaces and collinear magnetic systems. 
In this case, the spin-dependent charge and heat cur- 
rents are continuous across the interfaces and are scalar 



quantities. Past experiments show that whenever ox- 
ide layers are formed at interfaces the spin-dependent 
effects can be greatly enhancecpS. Furthermore, a non- 
coUinear system is required to describe important ap- 
plications of spi n-depe ndent transport such as the spin- 
torque-osdllatopSElHSS] ^j. gpin-transfer-torque magnetic 
memorySSMMl, Jn these cases the diffusion theory devel- 
oped here is not sufficient to describe the relevant pro- 
cesses. Instead, it should be described by a more general 
theory which includes spin-dependent tunneling and a 
3-dimensional spin vector, for example the magnetoelec- 
tronic circuit theorjW^'iSi Heikkila, Hatami and c owork- 
ers have previously developed such a theorjffSMl jn ex- 
plaining thermal spin-transfer torque and spin-dependent 
heat. 

In such a theory, the transport of electrons across inter- 
faces is described by a 4x4 conductance matrix G, which 
relates the total flux J = {Jc,Js,Q,Qs) at the inter- 
face to the variables at both sides u* = {V^,Vs,T^,Ts) 
(i=F,N) by J = G{u^-u^). 

Initial calculations on the conductance matrix for 
F/N interfaces have been carried out by Hatami and 
coworkers^ who in this framework calculated ther- 
mal spin-transfer-torque for various ferromagnet/non- 
magnetic interfaces. 

In the extended magnetoelectronic circuit theory they 
have introduced, all spin-dependent physics such as 
electrical or thermal spin injection is fully determined 
by the elements of the conductance matrix instead of 
the previously defined bulk spin-polarized parameters 
Pi,Ps,Pn,PQ of the ferromagnet. The physical origin 
in both theories is principally different; diffusion theory 
relates spin-dependent effects occurring at the interface 
to bulk ferromagnetic parameters (which in itself is only 
relevant very close to the interface) and the magnetoelec- 
tronic circuit theory describes the physics to occur due 
to quantum mechanical tunneling of spin states. How- 
ever, the effects they describe are the same. For exam- 
ple, thermal spin injection still describes the injection 
of spins in a non-magnetic metal due a heat flow over 
a ferromagnetic/non-magnetic interface. Dependent on 
the application, either the more simple diffusive theory 
developed here can be used or one needs to refer to the 
full circuit theory. 



VIII. SUMMARY 

We have developed a diffusive theory for spin- 
dependent charge and heat conduction which includes 
spin-orbit effects. Finite-element methods were used to 
model several experiments from literature where several 
parameters of this model were quantifled. Electrical spin 
injection, the spin-Hall angle of platinum and thermal 
spin injection were calculated. Also, new experiments 
were proposed which should demonstrate the spin-Peltier 
effect and a lower limit was given in an experiment which 
failed to demonstrate the magnetic heat valve. 
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